## [1] "explicated variable of regression : canopy_cover"
## [1] "for  Guinean_forest-savanna_regression_canopy_cover.RDS"           
## [2] "for  Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [3] "for  Sahelian_Acacia_savanna_regression_canopy_cover.RDS"          
## [4] "for  Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [5] "for  West_Sudanian_savanna_regression_canopy_cover.RDS"            
## [6] "for  Western_Congolian_forest-savanna_regression_canopy_cover.RDS" 
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "below, stancode for Guinean_forest-savanna_regression_canopy_cover.RDS"
## // generated with brms 2.20.4
## functions {
##   /* zero-inflated beta log-PDF of a single response
##    * Args:
##    *   y: the response value
##    *   mu: mean parameter of the beta distribution
##    *   phi: precision parameter of the beta distribution
##    *   zi: zero-inflation probability
##    * Returns:
##    *   a scalar to be added to the log posterior
##    */
##    real zero_inflated_beta_lpdf(real y, real mu, real phi, real zi) {
##      row_vector[2] shape = [mu * phi, (1 - mu) * phi];
##      if (y == 0) {
##        return bernoulli_lpmf(1 | zi);
##      } else {
##        return bernoulli_lpmf(0 | zi) +
##               beta_lpdf(y | shape[1], shape[2]);
##      }
##    }
##   /* zero-inflated beta log-PDF of a single response
##    * logit parameterization of the zero-inflation part
##    * Args:
##    *   y: the response value
##    *   mu: mean parameter of the beta distribution
##    *   phi: precision parameter of the beta distribution
##    *   zi: linear predictor for zero-inflation part
##    * Returns:
##    *   a scalar to be added to the log posterior
##    */
##    real zero_inflated_beta_logit_lpdf(real y, real mu, real phi, real zi) {
##      row_vector[2] shape = [mu * phi, (1 - mu) * phi];
##      if (y == 0) {
##        return bernoulli_logit_lpmf(1 | zi);
##      } else {
##        return bernoulli_logit_lpmf(0 | zi) +
##               beta_lpdf(y | shape[1], shape[2]);
##      }
##    }
##   // zero-inflated beta log-CCDF and log-CDF functions
##   real zero_inflated_beta_lccdf(real y, real mu, real phi, real zi) {
##     row_vector[2] shape = [mu * phi, (1 - mu) * phi];
##     return bernoulli_lpmf(0 | zi) + beta_lccdf(y | shape[1], shape[2]);
##   }
##   real zero_inflated_beta_lcdf(real y, real mu, real phi, real zi) {
##     return log1m_exp(zero_inflated_beta_lccdf(y | mu, phi, zi));
##   }
## }
## data {
##   int<lower=1> N;  // total number of observations
##   vector[N] Y;  // response variable
##   int<lower=1> K;  // number of population-level effects
##   matrix[N, K] X;  // population-level design matrix
##   int<lower=1> Kc;  // number of population-level effects after centering
##   int prior_only;  // should the likelihood be ignored?
## }
## transformed data {
##   matrix[N, Kc] Xc;  // centered version of X without an intercept
##   vector[Kc] means_X;  // column means of X before centering
##   for (i in 2:K) {
##     means_X[i - 1] = mean(X[, i]);
##     Xc[, i - 1] = X[, i] - means_X[i - 1];
##   }
## }
## parameters {
##   vector[Kc] b;  // regression coefficients
##   real Intercept;  // temporary intercept for centered predictors
##   real<lower=0> phi;  // precision parameter
##   real<lower=0,upper=1> zi;  // zero-inflation probability
## }
## transformed parameters {
##   real lprior = 0;  // prior contributions to the log posterior
##   lprior += student_t_lpdf(Intercept | 3, 0, 2.5);
##   lprior += gamma_lpdf(phi | 0.01, 0.01);
##   lprior += beta_lpdf(zi | 1, 1);
## }
## model {
##   // likelihood including constants
##   if (!prior_only) {
##     // initialize linear predictor term
##     vector[N] mu = rep_vector(0.0, N);
##     mu += Intercept + Xc * b;
##     mu = inv_logit(mu);
##     for (n in 1:N) {
##       target += zero_inflated_beta_lpdf(Y[n] | mu[n], phi, zi);
##     }
##   }
##   // priors including constants
##   target += lprior;
## }
## generated quantities {
##   // actual population-level intercept
##   real b_Intercept = Intercept - dot_product(means_X, b);
## }
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] " "
## [1] " "
## [1] "Guinean_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 1725) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -1.95      0.06    -2.06    -1.84 1.00     2329     2393
## fire_freq_std      -0.03      0.02    -0.07     0.02 1.00     2512     2410
## mean_precip_std     0.12      0.03     0.07     0.17 1.00     2246     2185
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi     5.28      0.19     4.93     5.68 1.00     2300     2519
## zi      0.11      0.01     0.10     0.13 1.00     2349     2292
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 243) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -1.77      0.36    -2.49    -1.06 1.00     2062     2093
## fire_freq_std      -0.00      0.06    -0.13     0.11 1.00     2320     2300
## mean_precip_std     0.10      0.17    -0.24     0.41 1.00     2021     2090
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi     6.72      0.62     5.57     7.98 1.00     2172     2174
## zi      0.06      0.02     0.03     0.09 1.00     2462     2175
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 5563) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -3.89      0.04    -3.96    -3.81 1.00     2460     2389
## fire_freq_std       0.18      0.03     0.12     0.23 1.00     2266     2215
## mean_precip_std     0.85      0.06     0.73     0.97 1.00     2353     2288
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi    27.54      1.24    25.14    29.96 1.00     2436     2332
## zi      0.68      0.01     0.67     0.69 1.00     2484     2090
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 47) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -1.68      2.43    -6.50     3.17 1.00     2565     2352
## fire_freq_std      -0.04      0.15    -0.34     0.25 1.00     2875     2352
## mean_precip_std    -0.20      1.18    -2.55     2.05 1.00     2539     2409
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi     4.47      1.51     2.01     7.78 1.00     2263     2328
## zi      0.51      0.07     0.38     0.64 1.00     2388     2116
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "West_Sudanian_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 3277) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -3.52      0.05    -3.61    -3.42 1.00     2378     2285
## fire_freq_std       0.08      0.01     0.05     0.11 1.00     2402     2369
## mean_precip_std     0.90      0.04     0.82     0.98 1.00     2507     2256
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi     9.97      0.32     9.36    10.63 1.00     2315     2251
## zi      0.23      0.01     0.21     0.24 1.00     2245     2288
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
##  Family: zero_inflated_beta 
##   Links: mu = logit; phi = identity; zi = identity 
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std 
##    Data: table_region (Number of observations: 259) 
##   Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
##          total post-warmup draws = 2400
## 
## Population-Level Effects: 
##                 Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept          -1.91      0.33    -2.54    -1.25 1.00     2474     2412
## fire_freq_std      -0.13      0.05    -0.23    -0.03 1.00     2465     2251
## mean_precip_std     0.20      0.21    -0.22     0.58 1.00     2423     2348
## 
## Family Specific Parameters: 
##     Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi     3.51      0.39     2.79     4.30 1.00     2551     2455
## zi      0.22      0.03     0.17     0.28 1.00     2147     2145
## 
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).

## [1] "########################################"
## [1] "Gamma regressions for  Guinean_forest-savanna_regression_rh98.RDS"           
## [2] "Gamma regressions for  Northern_Congolian_Forest-Savanna_regression_rh98.RDS"
## [3] "Gamma regressions for  Sahelian_Acacia_savanna_regression_rh98.RDS"          
## [4] "Gamma regressions for  Southern_Congolian_forest-savanna_regression_rh98.RDS"
## [5] "Gamma regressions for  West_Sudanian_savanna_regression_rh98.RDS"            
## [6] "Gamma regressions for  Western_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "Guinean_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 1725 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 8.17064
## [1] "sd(table_region$rh98)"
## [1] 4.228092

## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.061
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.058

## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.1
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.143

## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.231
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.174

## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.561
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.492

## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.023
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.201

## [1] "Northern_Congolian_Forest-Savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]  243 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 8.251706
## [1] "sd(table_region$rh98)"
## [1] 3.52155

## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.045
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.504

## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 7.762
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.474

## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.521
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.923

## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 7.937
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.17

## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.253
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.75

## [1] "Sahelian_Acacia_savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 5563 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 3.049537
## [1] "sd(table_region$rh98)"
## [1] 1.187586

## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.057
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.963

## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.037
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.972

## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.056
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 1.017

## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.057
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.999

## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.076
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 1.035

## [1] "Southern_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]   47 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 3.938866
## [1] "sd(table_region$rh98)"
## [1] 2.086246

## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 4.539
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 2.147

## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 4.376
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 2.142

## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 3.702
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 0.973

## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 3.446
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 1.511

## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 5.172
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 3.096

## [1] "West_Sudanian_savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 3277 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 5.728088
## [1] "sd(table_region$rh98)"
## [1] 3.15554

## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.838
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.069

## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.853
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.239

## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.647
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 2.967

## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.742
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.096

## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.674
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.021

## [1] "Western_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]  259 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$rh98)"
## [1] 6.137749
## [1] "sd(table_region$rh98)"
## [1] 4.607916

## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.332
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.828

## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.284
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 4.183

## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.204
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.881

## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.442
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.943

## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 5.531
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.084

## [1] "Beta regressions for  Guinean_forest-savanna_regression_canopy_cover.RDS"           
## [2] "Beta regressions for  Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [3] "Beta regressions for  Sahelian_Acacia_savanna_regression_canopy_cover.RDS"          
## [4] "Beta regressions for  Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [5] "Beta regressions for  West_Sudanian_savanna_regression_canopy_cover.RDS"            
## [6] "Beta regressions for  Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "Guinean_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 1725 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## Le chargement a nécessité le package : gtools
## 
## Attachement du package : 'gtools'
## Les objets suivants sont masqués depuis 'package:brms':
## 
##     ddirichlet, rdirichlet

## [1] "mean(table_region$canopy_cover)"
## [1] 0.1314814
## [1] "sd(table_region$canopy_cover)"
## [1] 0.14466

## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.162
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.15

## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.168
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.152

## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.167
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.158

## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.159
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.148

## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.153
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.144

## [1] "Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]  243 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$canopy_cover)"
## [1] 0.1649144
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1334132

## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.166
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.143

## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.199
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.158

## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.169
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.138

## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.158
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.125

## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.161
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.127

## [1] "Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 5563 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$canopy_cover)"
## [1] 0.005161202
## [1] "sd(table_region$canopy_cover)"
## [1] 0.02202149

## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.464
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.319

## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.478
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.317

## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.477
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.314

## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.472
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.317

## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.485
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.315

## [1] "Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]   47 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$canopy_cover)"
## [1] 0.04984602
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1063901

## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.393
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.256

## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.339
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.29

## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.435
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.283

## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.378
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.305

## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.358
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.297

## [1] "West_Sudanian_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 3277 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$canopy_cover)"
## [1] 0.05300694
## [1] "sd(table_region$canopy_cover)"
## [1] 0.08641946

## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.119
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.13

## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.117
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.126

## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.11
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.123

## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.115
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.127

## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.117
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.129

## [1] "Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1]  259 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"

## [1] "mean(table_region$canopy_cover)"
## [1] 0.1079882
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1649516

## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.168
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.19

## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.19
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.182

## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.187
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.185

## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.179
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.195

## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.187
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.192